        subroutine solvepseudo(Ainv,dudq,N_eq,numpar)

        real*8 Ainv(N_eq,N_eq)
        real*8 dudq(N_eq,numpar)
        real*8 sol(1000,100)


        do i=1,N_eq
            do j=1,numpar
                sol(i,j) = 0.
            enddo
        enddo        

        do k=1,numpar
            do i=1,N_eq
                do j=1,N_eq
                    sol(i,k) = sol(i,k) +  Ainv(i,j)*dudq(j,k)
                enddo
            enddo
        enddo

!         do k=1,numpar
!             r = 0.
!             do i=1,N_eq
!                 r = r + sol(i,k)**2
!             enddo
! !             write(*,*)'(x = A-1 * b) norma de x=  ',r
!         enddo

        do k=1,numpar
            do i=1,N_eq
                dudq(i,k)=sol(i,k)
            enddo
        enddo	

        return
        end
